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The more precise extraction for the CKM matrix element |Vc6| in the heavy quark 
effective field theory (HQEFT) of QCD is studied from both exclusive and inclusive 
semileptonic B decays. The values of relevant nonperturbative parameters up to order 
^ , l/ m Q are estimated consistently in HQEFT of QCD. Using the most recent experimental 

^s^j ■ data for B decay rates, \V c b\ is updated to be \V c b\ = 0.0395 ±0.0011exp ±0.0019^ n from 

B -v D*lv decay and \V cb \ = 0.0434 ± 0.0041 ex p ± 0.0020 th from B -> Dlu decay as well 
O ! as \V cb \ = 0.0394 ± 0.0010exp ± 0.0014 th from inclusive B -» X c lu decay. 

O: 

(S| ; PACS numbers: 12.15.Hh, 12.39.Hg, 13.20.He 

CN . Keywords: \V c b\, heavy quark effective field theory 

I. INTRODUCTION 

VO ■ 

■ The Cabibbo-Kobayashi-Maskawa (CKM) matrix element V c b describes the rich phenomena of flavor- 
changing transitions between the two heavy quarks b and c. Its precise extraction has become a very 
important issue in heavy flavor physics. Generally \V c b\ is extracted by studying either exclusive or 
inclusive semileptonic B decays. Since long distance contributions are involved in these decays, theoretical 
estimates of the nonperturbative parameters are of crucial importance for a precise extraction of \V c b\. 
Heavy quark symmetry (HQS) [1-3] and effective field theories of heavy quarks play a special role in such 

Qh! estimates. 

In this short note, we are going to provide a more precise extraction for |V^&| from both exclusive and 
| inclusive semileptonic B decays within the framework of HQEFT of QCD. This effective field theory can 

■ directly be derived from QCD [4] by integrating out the small components but with carefully treating 
the quark and antiquark components. Recently, it has been shown that the HQEFT can actually be 
regarded as a large component QCD [5]. Of particular, this effective field theory has been put forward 
and successfully applied to various hadron processes [6-16]. The resulting effective Lagrangian and the 
heavy quark expansion (HQE) in the HQEFT of QCD appear to be different from the usual heavy quark 
effective theory (HQET) [17,18] in which the quark and antiquark components were dealt with separately 
i.e., there is no quark-antiquark coupling terms considered in HQET. The application of the usual HQET 
in extracting \V c b\ has been discussed by several groups and summarized in ref. [19]. We shall further 
emphasize in this note the relevant new features of HQEFT with respect to the usual HQET . Our paper 
is organized as follows: \V c b\ extraction from exclusive and inclusive B decays will be discussed in section 
II and section III, respectively. The previously obtained values for \V c b\ will be updated by using the 
most recent experimental measurements. And a brief conclusion will be presented in section IV. 
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II. | Vcs | FROM EXCLUSIVE DECAYS 

The B — > D*(D)lv differential decay rates are 



with 
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^(1) =77^(1) =tm(1 + 6*), (2.3) 
Q{1) = nv[h + {l) mB ~ mD h-{l)] = W (l + (2-4) 



where the QCD radiative corrections to two loops give the short distance coefficients tja = 0.960 ± 0.007 
and r/v = 1.022 ±0.004 [20]. u> is the product of the four- velocities of the B and D*{D) mesons, u> = v-v'. 
The weak transition form factors /^(oj), h + {oj) and h-{oj) can be expanded in powers of 1/toq and 
represented by the heavy quark spin-flavor independent wave functions. Then based on Eqs.(2.1) and 
(2.2), \V c b\ can be precisely extracted as long as the form factors are reliably evaluated in some theoretical 
framework. 

The HQE of the transition form factors are studied in detail up to order 1 /to.q in HQEFT framework 
in Ref. [6]. When the contributions of operators containing two gluon field strength tensors are omitted, 
at the zero recoil point u = 1 the relevant form factors can be written as 

h Al = 1 + + 3« 2 ) - -L^ - n 2 )f - + 3F 2 



-2A gi - 6A^ 2 ) - —^-j- (Fi -F 2 - 2A gi + 2Ag 2 ) 



+ ^J- K , i F 1+ F 2 -2A gi -2A Q2) , (2.5) 

h+ = 1 + -^(- -) 2 [(«i + 3k 2 ) 2 - (F 1 + 3F 2 ) + 2A( gi + 3f> 2 )], (2.6) 

8A Z m& m c 

h- = 0, (2.7) 

where the parameters in the rhs. of Eqs.(2.5)-(2.7) are defined in HQEFT [6]. For simplicity, when 
the variable u) is not written explicitly, we refer to the zero recoil values of relevant functions, i.e., 
hAi = ^Ai(l)) K i = e tc The binding energy of a heavy meson M is defined as 

A = lim Am = lim (mM—mn). (2.8) 

It is seen from Eqs.(2.5) and (2.6) that automatically both the form factors h Al and h+ in HQEFT 
of QCD do not receive 1/mg order correction at the zero recoil point. Another point to be emphasized 
is Eq.(2.7). The vanished value of h- in HQEFT arises from the fact of partial cancellation between 
the 1/toq correction in the current expansion and the 1/toq correction coming from the insertion of 
the effective Lagrangian into the transition matrix elements [6]. Such a cancellation is not observed 
in HQET because in the latter framework the quark-antiquark couplings are not taken into account 
explicitly Though the so-called Luke's theorem in HQET protects the weak transition matrix elements 
from 1/toq order correction at zero recoil point, it does not protect h- from such correction [21]. As a 
result, in HQET framework, only the B — > D*lv decay rate at zero recoil is strictly protected against 
1/rriQ order correction while B — ► Dlv decay rate is not. This is the main reason besides the experimental 
considerations for the conclusion that the B — ► Dlv decay is not as favorable as B — > D*lv decay for 
\V c b\ extraction in HQET. 

In HQEFT, Eq.(2.7) advocates reliable extraction of |V^,| from both B — > D*lv and B — ► Dlv decays, 
because both rates of these decays do not receive 1/mg order correction, as can be seen from Eqs.(2.1)- 
(2.7). 

However, we note that the current world averages for | V^b|^"(l) and |T4i,|£(l) are [19] 

\V cb \T{l) = 0.0383 ± 0.0005 ± 0.0009, (2.9) 
\V cb \G{l) = 0.0413 ± 0.0029 ± 0.0027, (2.10) 

so the current experimental data for |14i>|£/(l) receive larger errors than those for |T4b|.F(l), which may 
lead to large experimental uncertainty for the result of \V c b\ extracted from B — > Dlv. 

Now the value of |T4f>| depends on our estimates of JF(1) and Q(l). First of all, it has been shown [6] 
that some of the nonperturbative parameters appearing in the rhs. of Eqs.(2.5) -(2.7) can be related to 
the heavy meson masses. Explicitly, we have 

A D(B) = A - (-J- - + 3« 2 ) - ^J-j^i + 3F 2 ) + 0{J-), (2.11) 

A^(^) = A - (J- - - «») - " F2) + 0i ^ } - (2 - 12) 

m c(b) Lm c (b) 4TO c(b) A m c(b) 



Thus some parameters can be determined from the heavy meson masses. A detailed study has been 
nresented in Ref. f6l. Here we would not renea.t the a.na.lvsis. Tt is easilv read from Ref. [fil that bv usins 



m b = 4.7GeV, m b - m c = 3.36GcV m b + A = 5.21GeV, 



(2.13) 



one obtains the following values for m, k 2 and Fi, F 2 : 

Ki « -0.615GeV 2 , k 2 « 0.056GeV 2 , 
F 1 « 0.917GeV 4 , F 2 w 0.004GcV 4 , (2.14) 

which agree well with the sum rule results: «i = -0.5±0.18GeV 2 , k 2 « 0.08GeV 2 and A = 0.53±0.08GeV 
[8]- 

Then up to order l/rriQ only two parameters g± and g 2 remain unknown. To have an estimation 
for them, we recall the definition of the nonperturbative parameters. They are defined by the matrix 
elements as follows [6]: 

< M^QlVr-^iip^Ql+^Mv >= - Kl (w)jTr[MTM} 

+ lTr[K aP {v,v')M l TP + l -a^M], 

< M^Q^T^—iip^i-iv ■ D)(ip ± )QM\M v >= - Qi (u)Itt[M'TM] 

iv ■ D A 

+ jTr[g al3 (v,v')M'TP + l -a^M}, 

< M^|Qit ) r-l^(z|) ± ) 2 -l^^) 2 Q(+)|M„ >= - Xl{L0 )^Tr[M'TM] 

+ ^Tr[ Xa p(v, v')M'TP + l -a aP M] + • • • , (2.15) 

where the ellipsis in the last equation represents the contributions of operators containing two gluon field 
strength tensors. The Lorentz tensor K a p is decomposed into scalar factors as K a p(v, v') — in 2 {uj)a a fj + 
K 3 (cj)(v' a ji3 — v'pj a ), and g a /3(v,v'), Xa[){v,v') are decomposed similarly. The functions F\ and F 2 are 
defined as 

Fi=xi+2Aei, (2.16) 
F 2 = X 2 + 2Ag 2 . (2.17) 

As will be further emphasized in the next section, the operator iv ■ D in the matrix elements was 
considered to give contribution of order the binding energy [6,7], i.e., iv ■ D ~ v ■ k ~ A , which has also 
been confirmed in the simple case by sum rule calculation [8]. Now if taking this simple replacement, we 
get from Eqs.(2.15) 

qi w —ki ■ A w 0.3GcV 3 , xi « «i « 0.4GcV 4 . (2.18) 

Note that our estimates in Eqs.(2.14) and (2.18) are consistent with the relation among them in Eq.(2.16). 

In Eq.(2.14) F 2 almost equals zero. Of course one reason for this may be the possible partial cancel- 
lation between \i an d 2Ag 2 in Eq.(2.17). But more importantly, we notice that \2, Q2 (and therefore 
F 2 ) are parameters characterizing the chromomagnetic type operators, and the contributions of these 
chromomagnetic operators are generally believed to be small. 

If using the same method of estimation as that for Q\ and \i, on e gets 

g 2 ~ — «2 • A w -0.03GeV 3 , (2.19) 
X 2 - k\ w 0.003GeV 4 . (2.20) 

It can be seen from the above equations that the two terms in Eq.(2.17) do have partial cancellation. 
However, Eqs.(2.19), (2.20) and (2.17) can not hold simultaneously. This indicates that Eqs.(2.19) and 
(2.20) are only rough estimates for the chromomagnetic type parameters since these parameters are very 
small in magnitude. Nevertheless, as g 2 must be small and can not influence the final result for \V cb \ 



1 This relation is in the sense of the effective contributions of operators within some matrix element. 



as significantly as g\ does, here we would first take Eq.(2.19) in extracting \V c b\ but leave more rigid 
determination of Q2 for future work. 

Fig.l presents 8* as a function of mj, + A at the fixed quark mass difference nib — m c = 3.36GeV. The 
curves are relatively flat in the region around rrib + A = 5.2GeV, which indicates a reliable extraction of 8* 
around this point. This value of nib + A is in agreement with that in Eq.(2.14). Varying mb and m& — m c 
in the ranges 4.6GeV < m b < 4.8GeV and 3.32GeV < m b - m c < 3.41GeV, and taking Sl w 0.3GeV 3 
and q 2 « -0.03GeV 3 , we get 

8* = 0.01 ± 0.04. (2.21) 
From Eqs.(2.9) and (2.21), we then obtain 

\V cb \ = 0.0395 ± O.OOllexp ± 0.001%. (2.22) 
For B — > Dlv decay, analogously, we obtain 

8 = -0.07 ±0.04 (2.23) 

and 

\V cb \ = 0.0434 ± 0.0041 ex p ± 0.0020 th . (2.24) 

Comparing Eqs.(2.22) and (2.24), we see that the results of \V c b\ extracted from B — ► D*{D)lv decays 
have similar theoretical errors. But large difference exists between the central values of our present results, 
Eqs.(2.22) and (2.24). However, due to the large experimental error for B — ► Dlv channel, the \V c b\ values 
extracted above are compatible with each other. A better determination of \V c b\ from B — > Dlv decay is 
expected when a more precise result for the quantity |T4f,|£(l) is obtained. 

Note that in Ref. [6] Q\ and Q2 are not estimated as above. There Q2 ~ O.lGeV 3 is assumed so that 
the central values for \V c b\ extracted from B — > D*{D)lv decays can be close to each other. Indeed, Q2 ~ 
O.lGeV 3 together with the experimental averages (2.9), (2.10) lead to \V cb \ = 0.0416±0.0011 ex p±0.0016 th 
from B -> D*lv decay and | V cb \ = 0.0417 ± 0.0040 ex p ± 0.0019 th from B -> Dlv decay. 



III. \V C b\ from inclusive decays 

Inclusive semileptonic B decays is the other alternative to determine \V c b\- In the usual HQET the light 
quark in a hadron is generally treated as a spectator, which does not affect the heavy hadron properties 
to a large extent. However, this treatment may be one possible reason for the failure of HQET in some 
applications. For example, the world average value for bottom hadron lifetime ratio T(Ab)/r(B () ) can not 
be explained well in the usual framework of HQET [22,23]. 

The dynamics of inclusive B decays is also analyzed in detail in Refs. [7,9] within the HQEFT frame- 
work. Instead of simply applying the equation of motion for infinitely heavy free quark, iv ■ DQ^ = 0, 
we treat the heavy quark in a hadron as a dressed particle, which means that the residual momentum 
k of the heavy quark within a hadron is considered to comprise contributions from the light degrees of 
freedom. This simple picture is adopted to conveniently take into account the effects of light degrees 
of freedom and the binding effects of heavy and light components of the hadron but not deal with the 
complex dynamics of hadronization directly. As the light degrees of freedom within the heavy hadron is 
relativistic, one has 

fc°~|k|. (3.1) 

Explicitly, the momentum of a heavy hadron H is represented as Pjj = itiqv + k + k' with k! being the 
momentum depending on the heavy flavor and suppressed by the inverse of the heavy quark mass. This 
directly leads to the relation between v ■ k and the binding energy, 

A = lim Ah — lim (run — mo) = v ■ k, (3.2) 

or 

(iv ■ D) = -= w A 



(3.3) 



To acquire a good convergence of HQE, we perform the expansion in terms of k — v(iv ■ D) (or say, 
equivalcntly, in terms of l/(rriQ + A)). Then the B — > X c lv decay rate is found to be [7,9] 



UB W = G l^J^ Vci(p,Pi^){Io(p,PiJ)+h(p,PiJ)^-h(p,PiJ)^}, (3.4) 



where fhb = nib + A, I , h an d h are functions of the mass square ratios p = ml/ml, p 2 = m 2 /m 2 
and pf = mf/rh 2 . Here the calculation is performed up to nonperturbative order l/rhg and perturbative 
order a 2 . The function rj c i characterizes QCD radiative corrections. Its two- loop results were obtained in 
Rcfs. [24,25]. rhb and m c can be determined from the meson masses via Eq.(2.11). k 2 is often extracted 
from the known B — B* mass splitting 

k 2 ~ \{m B , -m 2 B0 ) «0.06GcV 2 , (3.5) 
8 

which is consistent with Eq.(2.14) and the sum rule result [8]. 

There are several points to be mentioned for Eq.(3.4). Firstly, in deriving Eq.(3.4) the effects of light 
degrees of freedom are explicitly accounted for in the picture of a dressed heavy quark in a hadron. 
Secondly, it is seen that the next leading order contributions vanish in our HQE in terms of the inverse 
dressed heavy quark mass, 1/rhb- Furthermore, our HQE in terms of k — v{iv ■ D) (or 1/rhb) has a good 
convergence. It is found that the 1/ra 2 order contributions induce only —0.7 ~ 5% corrections to the 
total width T(Hb —* X c ev). Therefore we conclude that the higher order nonperturbative corrections can 
be safely neglected. Finally, now one needs only to treat the dressed quark mass rhb — "Tib + A instead of 
considering the uncertainties arising from the two quantities nib and A separately. Note that these are 
the features of HQE in HQEFT. They can not be observed in the HQE in the usual HQET, where one 
assumes (iv ■ D) to be zero or of higher order of 1/mb- In the latter case the next leading order corrections 
can be absent only when the HQE is performed in terms of l/mb, and the heavy quark mass mb and 
the binding energy A have to be treated separately. One result of this is that in HQET the theoretical 
prediction of the total decay width strongly depends on the value of bottom quark mass mj and may 
have larger uncertainties than in HQEFT. 

Using Eq.(3.4), \V c b\ can be extracted from experimental data for inclusive decay rates. Fig. 3 shows 
the obtained values of 5 in = T(B — > X c lv) / (\V c b\ 2 ) x 10 11 as a function of the energy scale p and the 
parameters m c , K\. It is seen that the extracted value of \V c b\ depends on the energy scale p weakly. So 
the main uncertainties come from m c and «i. The curves in Fig. 3(b) and Fig. 3(c) have minimal value 
points, around which the values of 5i n are favorable because they are less sensitive to the variation of 
m c and K\. Varying m c and K\ in the regions 1.45GeV <m c < 1.85GeV, -0.8GeV 2 < Kl < -0.4GeV 2 , 
these favorable values of 8i n change in the range: 

S m = 2.88±0.20ps-\ (3.6) 

where the central value is obtained at m c = 1.65GeV and K\ = — 0.6GcV 2 . 

Then the B° lifetime t(B°) = 1.540±0.014ps and the most recent CLEO data for B — > X c ev branching 
ratio Br(B -> X c ev) = (10.49 ± 0.17 ± 0.43)% [26] yield 

\V cb \ = 0.0394 ± 0.0010exp ± 0.001%. (3.7) 

\V c b\ extracted using different values of branching ratios is presented in Fig. 4. As mentioned above, the 
values at m c w 1.65GeV should be favorable because in this region the curves become less sensitive to 
the mass m c . Interestingly, as shown in Fig.3(c), when choosing m c = 1.65GeV we find that the curve of 
5 in (or the resultant |V^,|) as a function of — m reaches its minimal (or maximal) point at m « — 0.6GeV 2 . 
This value for K\ is again in good agreement with that we used in section II and with that obtained from 
sum rule calculation [8] . 

Note that the mass m c discussed here arises from the charm quark propagator in the HQE of the 
matrix elements. It should be the pole mass of charm quark. Thus it is not surprising that the value of 
m c obtained here is larger than the value for constituent mass quoted in section II. 



IV. CONCLUSION 



We have presented a more precise extraction of \V c b\ by studying the exclusive and inclusive semileptonic 
B decays up to the order of 1/toq in HQEFT of QCD. The nonperturbative parameters in HQEFT up 
to the same order are estimated consistentlv from various considerations. 



It has been shown that in HQEFT of QCD \V c b\ can be reliably extracted from exclusive decays 
B — > D*(D)lv with similar theoretical uncertainties, because neither of their differential decay rates 
receives 1/toq order corrections at zero recoil point. In studying inclusive B decays, we treat the heavy 
quark in a hadron as a dressed particle whose residual momentum comprises some effects from the light 
degrees of freedom. This enables us to consider the effects of light components in the hadron but still has 
a simple physical picture in application. In this framework, the B — > X c lv total decay rate is protected 
from l/m& order correction, and the 1/m^ order correction is very small. Furthermore, mj, and A only 
appear in the form of dressed quark mass: nib = nib + A, which must also reduce the uncertainties in our 
calculations as nib = mn [1 + 0(1 /nig)]. 

Using the most recent experimental data for the B-meson decay rates, we have arrived at the determi- 
nation for \V c b\ with 

|V c6 | = 0.0395 ± O.OOllcxp ± 0.0019 th from B D*lu, 
\V cb \ = 0.0434 ± 0.0041 CX p ± 0.0020 th from B -» Dlu, 
\V cb \ = 0.0394 ± 0.0010 e xp ± 0.0014 th from B X C W, 

where the result extracted from B — > Dlv decay receives a larger experimental uncertainty than that from 
B — ► D*lv decay but a similar theoretical uncertainty as the latter. The result obtained from B — ► D*lv 
decay agrees quite well with that from inclusive B — > X c lv decay. 
These results then give the average 

\V cb \ = 0.0402 ± 0.0014 cxp ± 0.0017 th 

Alternatively it can be represented as 

A = 0.83 ± 0.07 (4.2) 
in the Wolfenstein parameterization \V c b\ = A\ 2 with A = \V US \ = 0.22. 



(4.1) 
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FIG. 2. IVcbl extracted'frQj^ 1^ — ► D*/^ decay. (a),(c) are evaluated at mj — m c = 3.36GeV; (b) is evaluated 
at m b = 4.7GeV. 
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FIG. 3. Extracted 5 in = Fc(i? -> X c ei/)/(| \4 6 j 2 ) x 10 11 as a function of fi, m c and ki. (a): m = -0.6GeV 2 ; (b) 
and (c): n= 2.4GeV. 



FIG. 4. \V cb \ extracted from B — ► X c ev decay at fi = 2.4GeV and ki = -0.6GeV 2 



